Conditional Cooperativity in RAS Assembly Pathways on Nanodiscs and Altered GTPase Cycling.

Self-assemblies (i.e., nanoclusters) of the RAS GTPase on the membrane act as scaffolds that activate downstream RAF kinases and drive MAPK signaling for cell proliferation and tumorigenesis. However, the mechanistic details of nanoclustering remain largely unknown. Here, size-tunable nanodisc platforms and paramagnetic relaxation enhancement (PRE) analyses revealed the structural basis of the cooperative assembly processes of fully processed KRAS, mutated in a quarter of human cancers. The cooperativity is modulated by the mutation and nucleotide states of KRAS and the lipid composition of the membrane. Notably, the oncogenic mutants assemble in nonsequential pathways with two mutually cooperative 'α/α' and 'α/ß' interfaces, while α/α dimerization of wild-type KRAS promotes the secondary α/ß interaction sequentially. Mutation-based interface engineering was used to selectively trap the oligomeric intermediates of KRAS and probe their favorable interface interactions. Transiently exposed interfaces were available for the assembly. Real-time NMR demonstrated that higher-order oligomers retain higher numbers of active GTP-bound protomers in KRAS GTPase cycling. These data provide a deeper understanding of the nanocluster-enhanced signaling in response to the environment. Furthermore, our methodology is applicable to assemblies of many other membrane GTPases and lipid nanoparticle-based formulations of stable protein oligomers with enhanced cooperativity.

Identification of Potential Inhibitors Targeting GTPase-Kirsten RAt Sarcoma Virus (K-Ras) Driven Cancers via E-Pharmacophore-Based Virtual Screening and Drug Repurposing Approach.

BACKGROUND: Mutations in the K-Ras gene are among the most frequent genetic alterations in various cancers, and inhibiting RAS signaling has shown promising results in treating solid tumors. However, finding effective drugs that can bind to the RAS protein remains challenging. This drove us to explore new compounds that could inhibit tumor growth, particularly in cancers that harbor K-Ras mutations. METHODS: Our study used bioinformatic techniques such as E-pharmacophore virtual screening, molecular simulation, principal component analysis (PCA), extra precision (XP) docking, and ADMET analyses to identify potential inhibitors for K-Ras mutants G12C and G12D. RESULTS: In our study, we discovered that inhibitors such as afatinib, osimertinib, and hydroxychloroquine strongly inhibit the G12C mutant. Similarly, hydroxyzine, zuclopenthixol, fluphenazine, and doxapram were potent inhibitors for the G12D mutant. Notably, all six of these molecules exhibit a high binding affinity for the H95 cryptic groove present in the mutant structure. These molecules exhibited a unique affinity mechanism at the molecular level, which was further enhanced by hydrophobic interactions. Molecular simulations and PCA revealed the formation of stable complexes within switch regions I and II. This was particularly evident in three complexes: G12C-osimertinib, G12D-fluphenazine, and G12D-zuclopenthixol. Despite the dynamic nature of switches I and II in K-Ras, the interaction of inhibitors remained stable. According to QikProp results, the properties and descriptors of the selected molecules fell within an acceptable range compared to sotorasib. CONCLUSIONS: We have successfully identified potential inhibitors of the K-Ras protein, laying the groundwork for the development of targeted therapies for cancers driven by K-Ras mutations.

Structure-based prediction of Ras-effector binding affinities and design of "branchegetic" interface mutations.

Ras is a central cellular hub protein controlling multiple cell fates. How Ras interacts with a variety of potential effector proteins is relatively unexplored, with only some key effectors characterized in great detail. Here, we have used homology modeling based on X-ray and AlphaFold2 templates to build structural models for 54 Ras-effector complexes. These models were used to estimate binding affinities using a supervised learning regressor. Furthermore, we systematically introduced Ras "branch-pruning" (or branchegetic) mutations to identify 200 interface mutations that affect the binding energy with at least one of the model structures. The impacts of these branchegetic mutants were integrated into a mathematical model to assess the potential for rewiring interactions at the Ras hub on a systems level. These findings have provided a quantitative understanding of Ras-effector interfaces and their impact on systems properties of a key cellular hub.

Probing mutation-induced conformational transformation of the GTP/M-RAS complex through Gaussian accelerated molecular dynamics simulations.

Mutations highly affect the structural flexibility of two switch domains in M-RAS considered an important target of anticancer drug design. Gaussian accelerated molecular dynamics (GaMD) simulations were applied to probe the effect of mutations P40D, D41E, and P40D/D41E/L51R on the conformational transition of the switch domains from the GTP-bound M-RAS. The analyses of free energy landscapes (FELs) not only reveal that three mutations induce less energetic states than the wild-type (WT) M-RAS but also verify that the switch domains are extremely disordered. Principal component analysis (PCA) and dynamics analysis suggest that three mutations greatly affect collective motions and structural flexibility of the switch domains that mostly overlap with binding regions of M-RAS to its effectors, which in turn disturbs the activity of M-RAS. The analyses of the interaction network between GTP and M-RAS show that the high instability in hydrogen bonding interactions (HBIs) of GTP with residue 41 and Y42 in the switch domain I drives the disordered states of the switch domains. This work is expected to provide a molecular mechanism for deeply understanding the function of M-RAS and future drug design towards the treatment of cancers.

Dynamic regulation of RAS and RAS signaling.

RAS proteins regulate most aspects of cellular physiology. They are mutated in 30% of human cancers and 4% of developmental disorders termed Rasopathies. They cycle between active GTP-bound and inactive GDP-bound states. When active, they can interact with a wide range of effectors that control fundamental biochemical and biological processes. Emerging evidence suggests that RAS proteins are not simple on/off switches but sophisticated information processing devices that compute cell fate decisions by integrating external and internal cues. A critical component of this compute function is the dynamic regulation of RAS activation and downstream signaling that allows RAS to produce a rich and nuanced spectrum of biological outputs. We discuss recent findings how the dynamics of RAS and its downstream signaling is regulated. Starting from the structural and biochemical properties of wild-type and mutant RAS proteins and their activation cycle, we examine higher molecular assemblies, effector interactions and downstream signaling outputs, all under the aspect of dynamic regulation. We also consider how computational and mathematical modeling approaches contribute to analyze and understand the pleiotropic functions of RAS in health and disease.

Binding of active Ras and its mutants to the Ras binding domain of PI-3-kinase: A quantitative approach to KD measurements.

Ras family GTPases (H/K/N-Ras) modulate numerous effectors, including the lipid kinase PI3K (phosphatidylinositol-3-kinase) that generates growth signal lipid PIP3 (phosphatidylinositol-3,4,5-triphosphate). Active GTP-Ras binds PI3K with high affinity, thereby stimulating PIP3 production. We hypothesize the affinity of this binding interaction could be significantly increased or decreased by Ras mutations at PI3K contact positions, with clinical implications since some Ras mutations at PI3K contact positions are disease-linked. To enable tests of this hypothesis, we have developed an approach combining UV spectral deconvolution, HPLC, and microscale thermophoresis to quantify the KD for binding. The approach measures the total Ras concentration, the fraction of Ras in the active state, and the affinity of active Ras binding to its docking site on PI3K Ras binding domain (RBD) in solution. The approach is illustrated by KD measurements for the binding of active H-Ras and representative mutants, each loaded with GTP or GMPPNP, to PI3Kγ RBD. The findings demonstrate that quantitation of the Ras activation state increases the precision of KD measurements, while also revealing that Ras mutations can increase (Q25L), decrease (D38E, Y40C), or have no effect (G13R) on PI3K binding affinity. Significant Ras affinity changes are predicted to alter PI3K regulation and PIP3 growth signals.

Unorthodox regulation of the MglA Ras-like GTPase controlling polarity in Myxococcus xanthus.

Motile cells have developed a large array of molecular machineries to actively change their direction of movement in response to spatial cues from their environment. In this process, small GTPases act as molecular switches and work in tandem with regulators and sensors of their guanine nucleotide status (GAP, GEF, GDI and effectors) to dynamically polarize the cell and regulate its motility. In this review, we focus on Myxococcus xanthus as a model organism to elucidate the function of an atypical small Ras GTPase system in the control of directed cell motility. M. xanthus cells direct their motility by reversing their direction of movement through a mechanism involving the redirection of the motility apparatus to the opposite cell pole. The reversal frequency of moving M. xanthus cells is controlled by modular and interconnected protein networks linking the chemosensory-like frizzy (Frz) pathway - that transmits environmental signals - to the downstream Ras-like Mgl polarity control system - that comprises the Ras-like MglA GTPase protein and its regulators. Here, we discuss how variations in the GTPase interactome landscape underlie single-cell decisions and consequently, multicellular patterns.

Structure of the MRAS-SHOC2-PP1C phosphatase complex.

RAS-MAPK signalling is fundamental for cell proliferation and is altered in most human cancers1-3. However, our mechanistic understanding of how RAS signals through RAF is still incomplete. Although studies revealed snapshots for autoinhibited and active RAF-MEK1-14-3-3 complexes4, the intermediate steps that lead to RAF activation remain unclear. The MRAS-SHOC2-PP1C holophosphatase dephosphorylates RAF at serine 259, resulting in the partial displacement of 14-3-3 and RAF-RAS association3,5,6. MRAS, SHOC2 and PP1C are mutated in rasopathies-developmental syndromes caused by aberrant MAPK pathway activation6-14-and SHOC2 itself has emerged as potential target in receptor tyrosine kinase (RTK)-RAS-driven tumours15-18. Despite its importance, structural understanding of the SHOC2 holophosphatase is lacking. Here we determine, using X-ray crystallography, the structure of the MRAS-SHOC2-PP1C complex. SHOC2 bridges PP1C and MRAS through its concave surface and enables reciprocal interactions between all three subunits. Biophysical characterization indicates a cooperative assembly driven by the MRAS GTP-bound active state, an observation that is extendible to other RAS isoforms. Our findings support the concept of a RAS-driven and multi-molecular model for RAF activation in which individual RAS-GTP molecules recruit RAF-14-3-3 and SHOC2-PP1C to produce downstream pathway activation. Importantly, we find that rasopathy and cancer mutations reside at protein-protein interfaces within the holophosphatase, resulting in enhanced affinities and function. Collectively, our findings shed light on a fundamental mechanism of RAS biology and on mechanisms of clinically observed enhanced RAS-MAPK signalling, therefore providing the structural basis for therapeutic interventions.

Structure-function analysis of the SHOC2-MRAS-PP1C holophosphatase complex.

Receptor tyrosine kinase (RTK)-RAS signalling through the downstream mitogen-activated protein kinase (MAPK) cascade regulates cell proliferation and survival. The SHOC2-MRAS-PP1C holophosphatase complex functions as a key regulator of RTK-RAS signalling by removing an inhibitory phosphorylation event on the RAF family of proteins to potentiate MAPK signalling1. SHOC2 forms a ternary complex with MRAS and PP1C, and human germline gain-of-function mutations in this complex result in congenital RASopathy syndromes2-5. However, the structure and assembly of this complex are poorly understood. Here we use cryo-electron microscopy to resolve the structure of the SHOC2-MRAS-PP1C complex. We define the biophysical principles of holoenzyme interactions, elucidate the assembly order of the complex, and systematically interrogate the functional consequence of nearly all of the possible missense variants of SHOC2 through deep mutational scanning. We show that SHOC2 binds PP1C and MRAS through the concave surface of the leucine-rich repeat region and further engages PP1C through the N-terminal disordered region that contains a cryptic RVXF motif. Complex formation is initially mediated by interactions between SHOC2 and PP1C and is stabilized by the binding of GTP-loaded MRAS. These observations explain how mutant versions of SHOC2 in RASopathies and cancer stabilize the interactions of complex members to enhance holophosphatase activity. Together, this integrative structure-function model comprehensively defines key binding interactions within the SHOC2-MRAS-PP1C holophosphatase complex and will inform therapeutic development .

Structural basis for SHOC2 modulation of RAS signalling.

The RAS-RAF pathway is one of the most commonly dysregulated in human cancers1-3. Despite decades of study, understanding of the molecular mechanisms underlying dimerization and activation4 of the kinase RAF remains limited. Recent structures of inactive RAF monomer5 and active RAF dimer5-8 bound to 14-3-39,10 have revealed the mechanisms by which 14-3-3 stabilizes both RAF conformations via specific phosphoserine residues. Prior to RAF dimerization, the protein phosphatase 1 catalytic subunit (PP1C) must dephosphorylate the N-terminal phosphoserine (NTpS) of RAF11 to relieve inhibition by 14-3-3, although PP1C in isolation lacks intrinsic substrate selectivity. SHOC2 is as an essential scaffolding protein that engages both PP1C and RAS to dephosphorylate RAF NTpS11-13, but the structure of SHOC2 and the architecture of the presumptive SHOC2-PP1C-RAS complex remain unknown. Here we present a cryo-electron microscopy structure of the SHOC2-PP1C-MRAS complex to an overall resolution of 3 Å, revealing a tripartite molecular architecture in which a crescent-shaped SHOC2 acts as a cradle and brings together PP1C and MRAS. Our work demonstrates the GTP dependence of multiple RAS isoforms for complex formation, delineates the RAS-isoform preference for complex assembly, and uncovers how the SHOC2 scaffold and RAS collectively drive specificity of PP1C for RAF NTpS. Our data indicate that disease-relevant mutations affect complex assembly, reveal the simultaneous requirement of two RAS molecules for RAF activation, and establish rational avenues for discovery of new classes of inhibitors to target this pathway.

Therapeutic Targeting the Allosteric Cysteinome of RAS and Kinase Families.

Allosteric mechanisms are pervasive in nature, but human-designed allosteric perturbagens are rare. The history of KRASG12C inhibitor development suggests that covalent chemistry may be a key to expanding the armamentarium of allosteric inhibitors. In that effort, irreversible targeting of a cysteine converted a non-deal allosteric binding pocket and low affinity ligands into a tractable drugging strategy. Here we examine the feasibility of expanding this approach to other allosteric pockets of RAS and kinase family members, given that both protein families are regulators of vital cellular processes that are often dysregulated in cancer and other human diseases. Moreover, these heavily studied families are the subject of numerous drug development campaigns that have resulted, sometimes serendipitously, in the discovery of allosteric inhibitors. We consequently conducted a comprehensive search for cysteines, a commonly targeted amino acid for covalent drugs, using AlphaFold-generated structures of those families. This new analysis presents potential opportunities for allosteric targeting of validated and understudied drug targets, with an emphasis on cancer therapy.

Identification of functional substates of KRas during GTP hydrolysis with enhanced sampling simulations.

As the hub of major signaling pathways, Ras proteins are implicated in 19% of tumor-caused cancers due to perturbations in their conformational and/or catalytic properties. Despite numerous studies, the functions of the conformational substates for the most important isoform, KRas, remain elusive. In this work, we perform an extensive simulation analysis on the conformational landscape of KRas in its various chemical states during the GTP hydrolysis cycle: the reactant state KRasGTP·Mg2+, the intermediate state KRasGDP·Pi·Mg2+ and the product state KRasGDP·Mg2+. The results from enhanced sampling simulations reveal that State 1 of KRasGTP·Mg2+ has multiple stable substates in solution, one of which might account for interacting with GEFs. State 2 of KRasGTP·Mg2+ features two substates "Tyr32in" and "Tyr32out", which are poised to interact with effectors and GAPs, respectively. For the intermediate state KRasGDP·Pi·Mg2+, Gln61 and Pi are found to assume a broad set of conformations, which might account for the weak oncogenic effect of Gln61 mutations in KRas in contrast to the situation in HRas and NRas. Finally, the product state KRasGDP·Mg2+ has more than two stable substates in solution, pointing to a conformation-selection mechanism for complexation with GEFs. Based on these results, some specific inhibition strategies for targeting the binding sites of the high-energy substates of KRas during GTP hydrolysis are discussed.

Insights into the Cross Talk between Effector and Allosteric Lobes of KRAS from Methyl Conformational Dynamics.

KRAS is the most frequently mutated RAS protein in cancer patients, and it is estimated that about 20% of the cancer patients in the United States carried mutant RAS proteins. To accelerate therapeutic development, structures and dynamics of RAS proteins had been extensively studied by various biophysical techniques for decades. Although 31P NMR studies revealed population equilibrium of the two major states in the active GMPPNP-bound form, more complex conformational dynamics in RAS proteins and oncogenic mutants subtly modulate the interactions with their downstream effectors. We established a set of customized NMR relaxation dispersion techniques to efficiently and systematically examine the ms-µs conformational dynamics of RAS proteins. This method allowed us to observe varying synchronized motions that connect the effector and allosteric lobes in KRAS. We demonstrated the role of conformational dynamics of KRAS in controlling its interaction with the Ras-binding domain of the downstream effector RAF1, the first kinase in the MAPK pathway. This allows one to explain, as well as to predict, the altered binding affinities of various KRAS mutants, which was neither previously reported nor apparent from the structural perspective.

Chemical proteomic analysis of palmostatin beta-lactone analogs that affect N-Ras palmitoylation.

S-Palmitoylation is a reversible post-translational lipid modification that regulates protein trafficking and signaling. The enzymatic depalmitoylation of proteins is inhibited by the beta-lactones Palmostatin M and B, which have been found to target several serine hydrolases. In efforts to better understand the mechanism of action of Palmostatin M, we describe herein the synthesis, chemical proteomic analysis, and functional characterization of analogs of this compound. We identify Palmostatin M analogs that maintain inhibitory activity in N-Ras depalmitoylation assays while displaying complementary reactivity across the serine hydrolase class as measured by activity-based protein profiling. Active Palmostatin M analogs inhibit the recently characterized ABHD17 subfamily of depalmitoylating enzymes, while sparing other candidate depalmitoylases such as LYPLA1 and LYPLA2. These findings improve our understanding of the structure-activity relationship of Palmostatin M and refine the set of serine hydrolase targets relevant to the compound's effects on N-Ras palmitoylation dynamics.

Lipid Profiles of RAS Nanoclusters Regulate RAS Function.

The lipid-anchored RAS (Rat sarcoma) small GTPases (guanosine triphosphate hydrolases) are highly prevalent in human cancer. Traditional strategies of targeting the enzymatic activities of RAS have been shown to be difficult. Alternatively, RAS function and pathology are mostly restricted to nanoclusters on the plasma membrane (PM). Lipids are important structural components of these signaling platforms on the PM. However, how RAS nanoclusters selectively enrich distinct lipids in the PM, how different lipids contribute to RAS signaling and oncogenesis and whether the selective lipid sorting of RAS nanoclusters can be targeted have not been well-understood. Latest advances in quantitative super-resolution imaging and molecular dynamic simulations have allowed detailed characterization RAS/lipid interactions. In this review, we discuss the latest findings on the select lipid composition (with headgroup and acyl chain specificities) within RAS nanoclusters, the specific mechanisms for the select lipid sorting of RAS nanoclusters on the PM and how perturbing lipid compositions within RAS nanoclusters impacts RAS function and pathology. We also describe different strategies of manipulating lipid composition within RAS nanoclusters on the PM.

Blockade of mutant RAS oncogenic signaling with a special emphasis on KRAS.

RAS proteins (HRAS, KRAS, NRAS) participate in many physiological signal transduction processes related to cell growth, division, and survival. The RAS proteins are small (188/189 amino acid residues) and they function as GTPases. These proteins toggle between inactive and functional forms; the conversion of inactive RAS-GDP to active RAS-GTP as mediated by guanine nucleotide exchange factors (GEFs) turns the switch on and the intrinsic RAS-GTPase activity stimulated by the GTPase activating proteins (GAPs) turns the switch off. RAS is upstream to the RAS-RAF-MEK-ERK and the PI3-kinase-AKT signaling modules. Importantly, the overall incidence of RAS mutations in all cancers is about 19% and RAS mutants have been a pharmacological target for more than three decades. About 84% of all RAS mutations involve KRAS. Except for the GTP/GDP binding site, the RAS proteins lack other deep surface pockets thereby hindering efforts to identify high-affinity antagonists; thus, they have been considered to be undruggable. KRAS mutations frequently occur in lung, colorectal, and pancreatic cancers, the three most deadly cancers in the United States. Studies within the last decade demonstrated that the covalent modification of KRAS C12, which accounts for about 10% of all RAS mutations, led to the discovery of an adjacent pocket (called the switch II pocket) that accommodated a portion of the drug. This led to the development of sotorasib as a second-line treatment of KRASG12C-mutant non-small cell lung cancer. Considerable effort also has been expended to develop MAP kinase and PI3-kinase pathway inhibitors as indirect RAS antagonists.

Discovery of a dual Ras and ARF6 inhibitor from a GPCR endocytosis screen.

Internalization and intracellular trafficking of G protein-coupled receptors (GPCRs) play pivotal roles in cell responsiveness. Dysregulation in receptor trafficking can lead to aberrant signaling and cell behavior. Here, using an endosomal BRET-based assay in a high-throughput screen with the prototypical GPCR angiotensin II type 1 receptor (AT1R), we sought to identify receptor trafficking inhibitors from a library of ~115,000 small molecules. We identified a novel dual Ras and ARF6 inhibitor, which we named Rasarfin, that blocks agonist-mediated internalization of AT1R and other GPCRs. Rasarfin also potently inhibits agonist-induced ERK1/2 signaling by GPCRs, and MAPK and Akt signaling by EGFR, as well as prevents cancer cell proliferation. In silico modeling and in vitro studies reveal a unique binding modality of Rasarfin within the SOS-binding domain of Ras. Our findings unveil a class of dual small G protein inhibitors for receptor trafficking and signaling, useful for the inhibition of oncogenic cellular responses.

Spatiotemporal Imaging of Small GTPase Activity Using Conformational Sensors for GTPase Activity (COSGA).

Small GTPases cycle between active GTP bound and inactive GDP bound forms in live cells. They act as molecular switches and regulate diverse cellular processes at different times and locations in the cell. Spatiotemporal visualization of their activity provides important insights into dynamics of cellular signaling. Conformational sensors for GTPase activity (COSGAs) are based on the conserved GTPase fold and have been used as a versatile approach for imaging small GTPase activity in the cell. Conformational changes upon GDP/GTP binding can be visualized directly in solution, on beads, or in live cells using COSGA by fluorescence lifetime imaging microscopy (FLIM) technique. Herein, we describe the construction of COSGA for imaging K-Ras GTPase activity in live cells.

Using BioID to Characterize the RAS Interactome.

Identifying the proteins that associate with RAS oncoproteins has great potential, not only to elucidate how these mutant proteins are regulated and signal but also to identify potential therapeutic targets. Here we describe a detailed protocol to employ proximity labeling by the BioID methodology, which has the advantage of capturing weak or transient interactions, to identify in an unbiased manner those proteins within the immediate vicinity of oncogenic RAS proteins.

An Effective MM/GBSA Protocol for Absolute Binding Free Energy Calculations: A Case Study on SARS-CoV-2 Spike Protein and the Human ACE2 Receptor.

The binding free energy calculation of protein-ligand complexes is necessary for research into virus-host interactions and the relevant applications in drug discovery. However, many current computational methods of such calculations are either inefficient or inaccurate in practice. Utilizing implicit solvent models in the molecular mechanics generalized Born surface area (MM/GBSA) framework allows for efficient calculations without significant loss of accuracy. Here, GBNSR6, a new flavor of the generalized Born model, is employed in the MM/GBSA framework for measuring the binding affinity between SARS-CoV-2 spike protein and the human ACE2 receptor. A computational protocol is developed based on the widely studied Ras-Raf complex, which has similar binding free energy to SARS-CoV-2/ACE2. Two options for representing the dielectric boundary of the complexes are evaluated: one based on the standard Bondi radii and the other based on a newly developed set of atomic radii (OPT1), optimized specifically for protein-ligand binding. Predictions based on the two radii sets provide upper and lower bounds on the experimental references: -14.7(ΔGbindBondi)<-10.6(ΔGbindExp.)<-4.1(ΔGbindOPT1) kcal/mol. The consensus estimates of the two bounds show quantitative agreement with the experiment values. This work also presents a novel truncation method and computational strategies for efficient entropy calculations with normal mode analysis. Interestingly, it is observed that a significant decrease in the number of snapshots does not affect the accuracy of entropy calculation, while it does lower computation time appreciably. The proposed MM/GBSA protocol can be used to study the binding mechanism of new variants of SARS-CoV-2, as well as other relevant structures.